History and modes of star formation in the most active region of 
the Small Magellanic Cloud, NGC 346 1 
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ABSTRACT 

We discuss the star formation history of the SMC region NGC 346 based on Hubble Space 
Telescope images. The region contains both field stars and cluster members. Using a classical 
synthetic CMD procedure applied to the field around NGC 346 we find that there the star 
formation pace has been rising from a quite low rate 13 Gyr ago to w 1.4 x 10 -8 Mq yr~ 1 pc~ 2 
in the last 100 Myr. This value is significantly higher than in other star forming regions of 
the SMC. For NGC 346 itself, we compare theoretical and observed Color-Magnitude Diagrams 
(CMDs) of several stellar sub-clusters identified in the region, and we derive their basic evolution 
parameters. We find that NGC 346 experienced different star formation regimes, including a 
dominant and focused "high density mode" , with the sub-clusters hosting both pre-main sequence 
(PMS) and upper main sequence (UMS) stars, and a diffuse "low density mode" , as indicated by 
the presence of low-mass PMS sub-clusters. Quantitatively, the star formation in the oldest sub- 
clusters started about 6 Myr ago with remarkable synchronization, it continued at high rate (up 
to 2 x 10" 5 M Q yr 1 pc 2 ) for about 3 Myr and is now progressing at a lower rate. Interestingly, 
sub-clusters mainly composed by low mass PMS stars seem to experience now the first episode 
of star formation, following multi-seeded spatial patterns instead of resulting from a coherent 
trigger. Two speculative scenarios are put forth to explain the deficiency of UMS stars: the first 
invokes under-threshold conditions of the parent gas; the second speculates that the initial mass 
function (IMF) is a function of time, with the youngest sub-clusters not having had sufficient 
time to form more massive stars. 
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1. Introduction 

The Small Magellanic Cloud (SMC) is the clos- 
est star forming galaxy with a low metallicity {Z = 
0.004) typical of late- type dwarfs and most similar 
to that of primordial galaxies. For this reason an 
increasing number of studies are being devoted to 
its star formation history (SFH) and related pro- 
cesses. This research is part of a long-term project 
aimed at studying how the star formation started 
and propagated in the SMC, studying both young 
clusters and the field population. We concentrate 
here on the OB association NGC 346, the most 
active star forming region, where large numbers 
of pre-main sequence and massive stars are found. 
The inherent complexity of this star forming cloud 
is well recognized. This region provides an excel- 
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lent sample of newly formed stars spanning a wide 
range of masses, and bridging a wide range of tem- 
poral and spatial scales. 

The results presented in this paper are de- 
rived from observations acquired with the Ad- 
vanced Camera for Survey (ACS) on board of 
the HST. Several photometric investigations tack- 
led this complex population; all of them, using 
isochrone fitting, agreed that the star formation 
has taken place in a variety of sub-clusters at dif- 
ferent local conditions. Howe ver, different forma- 
tion scenarios are proposed: ISabbi et al.l (|2007l ) 
suggested a nearly coeval star formation in the 

(2000) 



Contursi et al 



cloud about 3 Myr ago. 
proposed a progressive star formation from the 
central cluster and pr opagating along the Bar. 
Gouliermis et al.l (|2008h argued that NGC 346 had 
been shaped by two delayed triggering events, the 
first one due to the central cluster and the second 
one to the massive progenitor of SNR B0057-724. 

The purpose of this study is to 1) re-examine 
the star formation of the individual star-forming 
sites to a much finer spatial scale (~ 1 pc), by 
means of the synthetic CMD approach, and 2) 
evaluate systematically whether and how the IMF 
and the star formation rate (SFR) are shaped and 
modified by the environment. 

2. Stellar content and its spatial distribu- 
tion in the region of NGC 346 

Fig. ED shows our HST image of NCG 346, 
acquired with three overl apping pointings of the 
ACS Wide Field Channel (|Nota et al.ll2006h . The 
data analys i s for this data set is presented in 
Sabbi et all (|2007l) (hereafter S07). Assuming an 
intrinsic distance modulus to the SMC of (m — 
M)o = 18.9, the field of view covers about 88 x 
88 pc 2 and contains different stellar populations, 
both clustered and diffuse. The clustered pop- 
ulation has two main co mponents: the interme - 
diate age cluster BS 90 (jBica fe Schmittl Il995l) . 
clearly recognizable as a roundish system at the 
top (north) of the image, and several clumps of 
stars ("sub-clusters" according to the definition 
by S07), spread over the field along a sort of 
umbrella shape, and identified by S07. The dif- 
fuse population is uniformly spread over the field. 
BS 90 is 4-5 Gyr old and most likely located at the 
foreground of the NGC 346 cluster (see S07 and 




Fig. 1. — ACS i mage for the ob served region 



around NGC 346 (|Nota et al.l l2006) 



Cignoni et ail2010l) . The sub-clusters are all part 



of the NGC 346 region, with ages f r om a few to 20 
Myr (S07, iGouliermis et al.l 120081 ICignoni et ail 



2010l ). while the diffuse population is presumably 



representative of the field SMC population in that 
region. The latter is then composed by SMC fore 
and background stars along the line-o f-sight, pos- 



sibly covering several kpc (see e.g. iGlatt et al 
120081 ). and spanning a wide range of ages, from 
Myr to ten Gyr. 

Our goal is to study the SFH of all these com- 
ponents by interpreting their observational CMDs 
with synthetic ones. BS 90 has already be en stud- 
ied by S07 (see also iRochau et al.l 120071) and no 



further analysis is presented here. We thus concen- 
trate on the other two components. The derivation 
of the SFH of the diffuse component is a stan- 
dard appli cation of the synthetic CMD method 
e.g., 



Tosi et al.l Il99li ICignoni et all 120061 



Cignoni fe Tosill2010h to a statistical y significant 



sample of field stars properly located in the re- 
gion far from the sub-clusters. On the other hand, 
the analysis of the clustered component is handled 
with a dedicated methodology for two main rea- 
sons: 1) the various sub-clusters appear to have 
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somewhat different ages and therefore we cannot 
treat them all together; 2) each sub-cluster con- 
tains few stars and therefore, when treated indi- 
vidually, has rather large associated statistical un- 
certainties. So, for the sub-clusters we have com- 
pared the observed CMDs with each other and 
with synthetic ones, taking into account the low 
number statistics and systematics. 

Following the class ification scheme proposed in 
Contursi et all (|2000l ). the NGC 346 region can be 



divided in a Spur, a filamentary low density struc- 
ture oriented to the North-East direction, and a 
fan-shaped feature (Bar), hosting the bulk of in- 
termediate and massive stars. 

The first question is whether two stellar tracers 
of the most recent activity, namely objects on the 
UMS and objects still on the PMS, are sharing the 
same "fine" spatial str ucture. A sim i lar an alysis 
has been performed bv lSchmeja et all (|2009h with 
the goal to derive the clustering behavior, while 
the purpose here is to identify sub-clusters with 
peculiar mass function and perform an accurate 
CMD analysis. 

The main advantage in using two different trac- 
ers stems from their different observational and in- 
trinsic uncertainties. With respect to UMS stars, 
PMS stars suffer less external contamination and 
low number statistics. On the other hand, UMS 
stars are less affected by incompleteness and the- 
oretical uncertainties. 

The regions of the CMD we associate with 
UMS and PMS stars are indicated in the top- 
left panel of Fig. [5] (with different colors in the 
electronic version). For UMS stars (black open 
diamonds) we chose all objects above the Turn- 
On (TOnj) of a 3 Myr isochrone and bluer than 
V -I = 0.2. With this definition, the UMS sam- 
ple is mostly composed by intermediate mass stars 
and by a few massive stars. For the PMS stars 
we considered two samples: objects redder than 
the 3 Myr isochrone with magnitude in the range 
22.5 < V < 25 (blue open circles, hereinafter 
PMS sample) and fainter than V = 25 (magenta 
open squares, low mass PMS sample, hereinafter, 
LPMS). It is easy to notice that while the PMS 
sample is younger than 3 Myr (or only slightly 



lr The TOn is the point in the color-magnitude diagram 
where the PMS st ars join the main sequence (see e.g. 
ICignoni et~aIll201Ch . 



older if a modest additional reddening is taken 
into account), the UMS sample can include MS 
stars as old as 600 Myr. Moreover, while there 
is no doubt that the PMS and LPMS samples are 
free from any contamination, the UMS sample can 
include a minor fraction of PMS stars starting to 
approach the MS. 

The top-right, bottom right and bottom left 
panels of Figure [2] show the location of the se- 
lected UMS, PMS and LPMS stars, respectively. 
These distributions provide a clue on the history 
of the region: half of the UMS stars are clumped 
into a few agglomerates (SC-1, SC-13 and SC-16, 
named according to S07), while the other half is 
more evenly distributed, and are probably mem- 
bers of the SMC field (foreground and background 
stars). One of the most striking aspects of this dis- 
tribution is the absence of filamentary structures 
and a rather round appearance of the sub-clusters. 

In contrast, the PMS stars are found almost 
exclusively either clumped or irregularly arranged 
along filaments. More filamentary than clumpy is 
also the distribution of the LPMS sample. In par- 
ticular, we note that the Spur region, composed 
by a few aggregates in the PMS map, becomes a 
sort of bridge extending for tens of pc in the LPMS 
map. A word of caution is however necessary for 
the frequent holes in their star distribution: like 
that clearly visible in the center of SC-1, they are 
likely due to incompleteness effects. Incomplete- 
ness is caused by crowding, which reaches a max- 
imum in the central region SC-1, and is much less 
severe at the outskirts of the region. 

3. Sub-cluster properties 

Looking at the maps of Figure [U there is a fur- 
ther intriguing aspect of the star spatial distribu- 
tion: for reasons that will be discussed later, not 
all the sub-clusters visible in the UMS map are de- 
tected in the PMS map, and vice versa. In Figure 
[H stellar sub-clusters which host both species are 
indicated with red arrows, while blue and green 
arrows indicate sub-clusters which host only PMS 
or UMS stars, respectively. We find that the PMS 
sub-clusters reside in the North-East part of the 
Spur (NE-Spur) (SC-14, SC-15 and SC-15b) and 
in the North- West side of the Bar (NW-Bar) (SC- 

4, SC-5 and SC-6). On the other hand, UMS-PMS 
sub-clusters are found in the central region (SC- 
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Fig. 2. — Top-left panel: CMD for the region NGC 346. Black open diamonds, blue open circles and 
magenta open squares indicate the areas referred to the UMS, P MS and LPMS star s, respectively. The 
line corresponds to a 3 Myr isochrone including the PMS phase ([Cignoni et al. The other panels 



(from top-right clockwise) show the corresponding spatial distributions. Red arrows indicate sub-clusters 
hosting UMS and PMS stars; blue arrows refer to PMS dominated sub-clusters; and the green arrow refers 
to the only UMS sub-cluster. The number labeling each sub-cluster follows the nomenclature provided by 



Sabbi et al.l (|2007l ). where the whole catalogue with the star coordinates can also be found. In the top-right 



panel we shows also the RA and DEC coordinates of the plot extremes and of its center (X=Y=4000). 



1), in the southern cluster of the Spur (SC-13) and 
in the South-East side of the Bar (SE-Bar) (SC-9, 
SC-10 and SC-11). The only sub-cluster to host 
exclusively UMS stars is SC-16. In fact, SC-16 is 



older than 10 Myr (see e.g. S07 and lCignoni et al 



2010), and its PMS members are too blue to be in- 
cluded in our PMS selection. The situation is rad- 
ically different for SC-14, SC-15 and SC-15b. De- 
spite the few Myr age of these sub-clusters, there 
are few, if any, UMS stars: where is the MS coun- 
terpart of their observed PMS stars? 

To try and understand the sub-cluster similari- 
ties and differences, we have divided them in three 
categories: Group I sub-clusters are those with 



minimum PMS/UMS ratio, Group II sub-clusters 
those with highest PMS/UMS ratio, while we have 
labelled as Group III the sub-clusters where the ra- 
tio was either intermediate or difficult to estimate. 
Fig. |3] shows the clear spatial separation of the 
three different Groups. Figure 2] shows the CMD 
for typical Group I sub-clusters (SC-1 and SC-13) 
and for Group II sub-clusters (SC-14 and SC-15) 
as a function of the distance from the sub-cluster 
center. Different rows refer to stars selected from 
four equal area annuli. 

Excluding the very central region of SC-1, af- 
fected by incompleteness, these CMDs trigger two 
interesting questions: 
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Fig. 3. — The three boxes superimposed to the NGC346 map identify sub-clusters with similar proper- 
ties. Solid circles indicate the sub-clusters whose CMDs are discussed in the text. The dashed large circle 
represents the region of field that has been used to estimate the contamination (see text). 



1) Lack of UMS stars: As suggested by the spa- 
tial maps, intermediate and massive MS stars seem 
deficient in sub-clusters like SC-14. Massive stars 
are always rare and sub-clusters like SC-14 are tiny 
pockets of stars; hence, stochastic fluctuations due 
to small number statistics in the poorly populated 
high mass end of the IMF need to be carefully ac- 
counted for. Is this sufficient to conclude that the 
star formation in SC-14/15 is a simple scaled-down 
version of that of SC-1? 

2) Redness of PMS stars: The red tail of the 
PMS distribution extends to redder colors going 
from SC-1 to SC-15 along the spur. Is their 
redness intrinsic (age dispersion or circumstcllar 
disks) or caused by obscuring material along the 
line of sight? 

From the analysis of the sub-clusters CMD we 
find that the Group I members show minimum 
PMS redness, minimum UMS spread and, by def- 
inition, minimum PMS/UMS ratio. The Group II 
members show the maximum PMS redness and, 
by definition, the highest PMS/UMS ratio. The 
Group III members resemble Group I members in 
relation to the PMS redness and to Group II in 
relation to the PMS/UMS ratio, while the large 
color spread among UMS stars is a unique prop- 
erty of this class of sub-clusters. 

The next section is dedicated to a quantita- 
tive analysis of these groups. In order to re- 



produce the observed CMD morphology and star 
counts of the sub-clusters shown in Figure [2] we 
applied the synthetic CMD approach, combining 
the Z = 0.004 Pisa stella r models for PMS stars 
(see ICignoni et al. 2009) with the same metal- 



licity Padova stellar models ([Fagotto et al.l 119941 ) 
for later evolutionary phases. The explored mass 
range is between 0.45 and 120 M©. To produce 
realistic simulations we have incorporated photo- 
metric errors and incompleteness corrections as 
derived from extensive artificial tests on the real 
images (S07). Different parameters, such as ex- 
tinction, IMF, binarity, and star formation rate 
are varied until the stellar densities and distribu- 
tions well match the observed CMDs. 

To better interpret the morphology of the sub- 
cluster sequences we chose sub-cluster radii as the 
best compromise between the need to bypass the 
crowding problems, typical of the very inner cen- 
tral regions, and the need to minimize the effect 
of field star contamination, while still containing a 
reasonable number of sub-cluster stars. Concern- 
ing SC-1, we find that the best strategy is to focus 
on an annulus around the center from 50 to 282 
pixels (see Figure [3]). Artificial star tests indicate 
that in this area stellar detections are more than 
90% complete down to V ~ 24. For all minor and 
satellite sub-clusters we find that a circle of radius 
141 pixels (see Figure [3]) allows to include most of 



5 



SC-1 SC-13 SC-14 SC-15 



V 19 



V 19 



V 19 



I I I I I I I I I I 1 I 1 I 1 

— • 
• 

% 

\ 


i I i I i I i I i I i I i I i 
• 

• 
• 

4 • •• 
• 

V*-. 
• •• • V 


I | I | I | I | I | I | I | I 
• 

• 

i - . 


I | I | I | I 

• 
• 

* 

• 

= « 


i i i i i i i 
• 


1' 1 ' 1 1 1 1 1 1 1 ' 1 ' 1 ' 

• 

• % 

• • • 


• 
• 

m 
• 

• 

• • 

-.hhi- 1 .-r.'i . i 


; ' 1 ' 1 1 1 1 1 > 1" 1 i 1 i 

• 

• 

• • a- • 

* 


• 


• 

• 

'* [ i ll 1 


i < i < i < i < i ' i ' i ' 

_ • 

— • 

- • 

: .**$ 

-- 


• 

• 

• 


• • 
• 
• 

• - • 

; rV . V i -i i i 


; i I i I ' I 



• • 


- • 


• 

-.I.T,T,f*T,T«T, 


• 

v • 

-,T,T,T,W,*T,«J.t, 





012 012 012 012 



V-l V-l V-l V-l 



Fig. 4. — From left to right, the four columns show CMDs for SC-1, SC-13, SC-14 and SC-15 respectively. 
From top to bottom, stars are shown that belong to progressively more external annuli of equal area centered 
on the highest density peak. The radius of the inner circle is about 71 pixels (rj 1 pc) for all sub-clusters. 



the members with a modest field contamination. 

4. Group I 
4.1. SC-1 

SC-1 is the central and most populous sub- 
cluster in NGC 346. We have simulated synthetic 
CMDs based on different combinations of evolu- 
tionary (SFH, IMF) and environmental parame- 
ters (reddening, distance modulus, fraction of bi- 
nary stars). The best combination of parameters 
is assessed by examining the star counts in strate- 
gic boxes along the main sequence (see labels MSI 
to MS4 in Figure top panel) and in the pre- 
main sequence (PMS). Given the youth of this sub- 



cluster, likely younger than 10 Myr, regions MSI 
and MS2 are particularly suitable to constrain the 
IMF. MS3 and MS4 convey information both on 
the IMF and on the star formation history. The 
PMS region informs mainly on the IMF (given the 
uncertainties on the PMS models, the PMS box 
is used only as, a posteriori, consistency check). 
From a numerical point of view, a grid search rou- 
tine is used to effectively determine the combina- 
tion of parameters minimizing residuals in these 
regions. To limit the parameter space we also 
considered two additional morphological features: 
the magnitude V dispersion of the PMS stars at 
the TOn (marked with arrows in the top panel 
of Figure [5]) , and the number of MS stars in the 
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range of magnitude 21.3 and 21.8 (approximately 
the 6 Myr TOn and the BS 90 Turn-Off), which 
is a strong indicator of contamination and/or sub- 
cluster members older than about 6 Myr. Our 
solution to handle this contamination was to use 
as control fields hundred regions of equal area lo- 
cated in the south-west part of the image (see the 
large dashed circle in Figure EJ). 

As a first result, a Salpeter IMF (jSalpeter 



19551 ). a fraction of binaries of 30%, a foreground 
reddening E(B — V) = 0.08 and a distance mod- 
ulus (m — M) = 18.9 provide a good agreement 
with the data. Concerning the SFH, we find that 
the stellar production in SC-1 started energeti- 
cally between 5 and 6 Myr ago, was strong for 
about 3 Myr (period between 3 and 6 Myr ago), 
then it dropped, probably victim of feedback from 
the massive stars of the first generation which 
quenched the subsequent formation. Before the 
onset at 6 Myr, the observed counts in the CMD 
window between the 6 Myr TOn and the BS 90 
Turn-Off suggest a null or negligible activity: out 
of 9 objects in the magnitude interval 21.3 — 21.8 
and V—I < 0.45 at least three can be attributed to 
the field with a confidence level better than 95%. 
The residual stars are compatible with a star form- 
ing activity of at most 0.3 x 10~ 5 Mq yr _1 pc~ 2 in 
the period between 8 Myr and 6 Myr ago. During 
the active phase this region experienced a pealo 
of about 2 x 10^ 5 M Q yr 1 pc 2 between 4 and 5 
Myr ago, a val ue much higher than in the center of 
NGC 602 Csee lCignoni et alJl2009h . another active 
star forming region in the SMC. After this onset, 
the MS boxes strongly exclude that the star for- 
mation was constant or increasing to the present 
day: according to our best model, only 20% of 
the total mass of young stars in the explored re- 
gion is produced in the last 3 Myr. Moreover, the 
smoothness of the UMS is not suggestive of any 
recent short burst. 

Figure [5] shows a comparison between our best 
synthetic CMD (bottom) with the observational 
one (top). For better visualization, the model has 
been artificially "contaminated" with field stars 
taken from an appropriate region of equal area. 
In broad terms, the UMS and the TOn regions 



2 Actually, this is a lower limit, which does not take into 
account stars below 0.45 Mq. When extrapolated using a 
Salpeter slope down to 0.1 Mq, this translates in an upper 
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Fig. 5. — Comparison of observational (top panel) 
and best fit synthetic CMD (bottom panel) for the 
sub-cluster SC-1. Also shown are the regions used 
for characteristic star counts (see text). There 
are four such regions along the MS and one for 
the PMS phase. The theoretical CMDs are cal- 
culated with a 30% of binaries and Salpeter IMF. 
In the field shown in the figure, the observational 
and theoretical CMDs contain the same number 
of stars. 



are well reproduced, although the synthetic main 
sequence appears more fuzzy. Moreover, two dis- 
crepancies are noted: 1) our models do not pro- 
vide satisfactory results about the dispersion of 
PMS stars, predicting smaller spreads and bluer 



color s than observed (see also iPozzo et al.l 12003 



limit of 3.6 X 10 -5 M, 



®yr x pc 



and Mavne et al. 20071 ): 2) no combination of pa- 
rameters reproduces the observed number of PMS 
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stars. Once the MS boxes are matched, our best 
model underestimates it systematically by about 
40%. Although part of this difference may arise 
from our lack of PMS models less massive than 
0.45 Mq (potentially present in the PMS box for 
ages younger than few hundreds of Kyr), a com- 
plete solution of the discrepancy seems to require 
other physical reasons. It is intriguing to note that 
such discrepancy could be easily resolved by con- 
sidering an additional episode of star formation 
younger than 2 Myr, accounting for about 500 Mq, 
where the MS phase is artificially suppressed and, 
therefore, not visible in the CMD. As an alter- 
native, the IMF could be steeper than Salpeter's 
for masses below 2 M© , thereby creating a larger 
population of low mass stars. 

In the following we adopt SC-1 as reference sub- 
cluster to compare with minor sub-clusters. 

4.2. SC-13 

The SC-13 sub-cluster is less dense than SC- 
1: excluding MS stars fainter than V — 22, it 
accounts for about 7 stars per pc 2 against about 
11 stars per pc 2 found in SC-1. Moreover, SC-13 
is less affected by incompleteness. Despite these 
differences, when the CMD of SC-1 is normal- 
ized to the number of SC-13 stars brighter than 
V = 22 (inside a radius of 141 pixels), the CMD 
morphologies appear very similar (see Figure [5]) ■ 
First, the color dispersion along the upper main 
sequence and the TOn magnitude are indistin- 
guishable from one another, suggesting an equiv- 
alent star formation duration and onset. Second, 
the PMS region displays an identical spread and 
the ratio PMS /UMS is fairly simila r (about 2, see 
Table [J) . As already envisaged in ICignoni et al 



( 2010h . the TOn gets brighter (up to 1 mag) when 
stars from the inner region of SC-13 (see the top 
row, second column of Figured]) are selected. Al- 
though numbers are too small to allow statistically 
significant conclusions, these findings suggest that: 
1) the star formation in SC-13 was triggered at 
the same time of SC-1; 2) the inner region of SC- 
13 may have experienced a secondary star forming 
episode about 3 Myr ago. In terms of rate density, 
we estimate that SC-13 astrated (at peak) about 
1.3 x KT 5 M Q yr- l pc- 2 . 
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Fig. 6. — Comparison between SC-1 (filled red cir- 
cles) and SC-13 (open blue diamonds). The CMD 
for SC-1 is normalized to the number of SC-13 
stars brighter than V = 22. 

Table 1: UMS, PMS, LPMS star counts and 
PMS/UMS ratio (P/U) for the labeled sub- 
clusters. All counts are measured inside a radius 
of 141 pixels, apart SC-1 that is measured in an 
annulus around the center between 50 and 282 pix- 
els. 







UMS 


PMS 


LPMS 


P/U 




SC-1 


94 


195 


79 


2.1 


Group I 


SC-13 


20 


35 


13 


1.8 




SC-9 


10 


19 


7 






SC-10 


7 


19 


1 


2.6 




SC-11 


7 


24 


18 






SC-14 


5 


23 


19 




Group II 


SC-15 





24 


24 


10.8 




SC-15b 





7 


11 






SC-4 


2 


22 


19 




Group III 


SC-5 


2 


29 


11 


10.4 




SC-6 


3 


22 


18 





4.3. SC-9, SC-10, SC-11 

In Figure [7] we present (blue open diamonds) 
the combination of the CMDs of the sub-clusters 
SC-9, SC-10 and SC-11, all located in the SE-Bar, 
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Fig. 7.— Composite CMD for SC-9, SC-10 and 
SC-11 (open blue diamonds) superimposed to the 
CMD for SC-1 (filled red circles) normalized to 
have the same number of stars brighter than V = 
22. 

superimposed to SC-1 (conveniently normalized 
to the same number of stars brighter than V = 
22). The similarity of morphologies is remark- 
able: UMS, intermediate mass PMS and PMS 
spread overlap. In terms of star counts, the ratio 
PMS/UMS (about 2.6, see Table [J) is consistent 
within the expected Poisson fluctuations. 

When normalized to the same area, we estimate 
the maximum star formation rate density in the 
SE-Bar to be 0.4 x lfr 5 M Q yr' 1 P c~ 2 . 

5. Group II 

All located in the NE-Spur, the members of this 
group are embedded in dust and nebulosities. Un- 
like Group I, these sub-clusters are mainly com- 
posed of low-mass stars with few, if any, UMS 
stars. This condition strongly hinders any attempt 
to obtain reliable ages using the TOn. Figure |8] 
shows a direct comparison between the compos- 
ite CMD of SC-14, 15 and 15b with SC-1 (re- 
sampled to have the same number of bright stars 
with V < 22). Two differences emerge: 1) the 
Group II CMD shows a ratio PMS/UMS stars of 
about 11 (see Table [1]), at least a factor of three 
larger than in SC-1, and a lack of intermediate 
mass PMS stars, which are clearly present in the 



Fig. 8. — Comparison between Group II sub- 
clusters SC-14, 15, 15b (open blue diamonds) and 
SC-1 (filled red circles) re-sampled to have the 
same number of stars brighter than V = 22. 

CMD of SC-1 at V = 20 and in the color range 
0.1 < V - I < 1; 2) In the Group II CMD the 
LPMS stars are much redder than in SC-1. 

Concerning point 1), we notice however that the 
presence of intermediate mass and massive young 
stellar objects (YSO) has been discovered with 
Spitzer observations by ISimon et al. (2007). sug- 
gesting that a fraction of UMS stars may exist but 
be still invisible in optical wavelengths. We show 
in Figure IHlthe distribution of Spitzer YSOs in our 
field of view. 

Point 2) is not caused by either dust between us 
and the sub-cluster or diffuse dust within the sub- 
cluster (otherwise also the lower main sequence 
would exhibit the color excess) but rather by red- 
dening material intimately related to the individ- 
ual PMS stars. It is also noteworthy that the UMS 
is not reddened: either these objects belong to 
the field or the reddening material affects only the 
PMS phase (as expected for circumstellar mate- 
rial) . 

Among the Group II sub-clusters, SC-15 (Fig- 
ure [TU]) is the one displaying more significant dif- 
ferences with respect to SC-1. The upper main 
sequence is definitely underpopulated, with just a 
couple of stars at best. In addition, the LPMS 
stars in SC-15 are at least 0.5 mag redder than 
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Fig. 9. — Spatial distribution of the Spitzer sample 
of intermediate an d high mass YSOs (data from 
ISimon et al.l l2007l) overlaid on the map of PMS 
(black dots) and UMS stars (large black dots) from 
Figure ^ . Red circles and green diamonds stand 
for YSOs with estimated mass between 4.5 and 8 
Mq and larger than 8 M Q , respectively. 



comparable objects in SC-1, while the main se- 
quence stars fainter than V = 22 are only slightly 
redder (0.1 mag). Once again, this would support 
the idea that reddening material is differentially 
distributed among PMS stars. 

In an attempt to reconcile the PMS samples, 
we shifted the SC-15 stars along the reddening 
vector until the color distribution of PMS in the 
two sub-clusters CMD overlapped. We derived 
the optimum shift when the two brighter stars in 
the SC-15 PMS sample reached the main sequence 
color. Figure [TO]- (bottom panel) shows the CMD 
obtained with this methodology. It is clear that 
the resulting correction, E(V — I) ~ 0.6, is still 
largely insufficient to reconcile the bulk of PMS 
stars. Only a differential reddening may overcome 
the remaining discrepancy. However, the redden- 
ing vector (see Figure [101 bottom panel) is al- 
most parallel to the PMS and a full (differential) 
correction would produce too many bright PMS 
stars. In other words, although a strong extinction 
by circumstellar envelopes would be very likely in 
such young objects, it would necessarily further 
increase the ratio of PMS stars over UMS stars. 



13 
15 
17 

V 19 

21 
23 
25 

13 
15 
17 

V 19 

21 
23 
25 



• SC-1 
o SC-15 



0.0 



0.5 



1.0 1.5 
V-l 



2.0 



V • 

• oo v * 40 
• «*o4°.« % •» 



E(V-I) 



Km 



0.0 



0.5 



1.0 1.5 
V-l 



2.0 



Fig. 10. — Top panel: comparison between SC-15 
(open blue diamonds) and SC-1 (filled red circles) 
normalized to the same area. Bottom panel: SC- 
15 after reddening correction applied to PMS stars 
only. 



6. Group III 

All members of this group are located in the 
NW-Bar. Figure [TTJ shows the combination of 
CMDs for the sub-clusters SC-4, SC-5 and SC-6 
superimposed to the CMD for SC-1. The stel- 
lar population in this group shows CMD features 
which are somehow intermediate between Group 
I and Group II. The color spread in the LPMS is 
fully consistent with the same sequence in SC-1. 
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Fig. 11. — Comparison between the composite 
CMD for Group III sub-clusters SC-4, SC-5 and 
SC-6 superimposed to the CMD for SC-1 normal- 
ized to have the same number of stars brighter 
than V = 22. 
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Fig. 12. — Spatial position of the selected field 
stars (red dots), PMS stars younger than 3 Myr 
(blue dots) and main sequence st ars with 23 < 
V < 24 (black dots; see Fig. 3 in Cignoni et al.l 
l201dt ). respectively. 



On the other hand, the ratio PMS/UMS of about 
10.4 (see Table [IJ is close to the value found in 
Group II. Moreover, while in SC-1 the UMS is 
sharply defined and an extended sequence of in- 
termediate mass PMS stars is observed, all bright 
(V < 23) stars in the Group III CMD seem not 
to lie along any r ecognizable s e quenc e. Also for 
these sub-clusters ISimon et al.l ([20071 ) report the 
detection of massive embedded YSOs (see Fig. [9]). 

Finally, it is worth noticing that the field con- 
tamination, represented by MS stars fainter than 
V = 22, is stronger in the NW-Bar and in the 
NE-Spur than in the SE-Bar because the former 
regions are closer to the cluster BS 90. 

7. Field Star Formation History 

In addition to a very young population, clearly 
reflected by clumps or filaments mainly studded 
by PMS stars, Figure [2] reveals the signatures of 
an evenly distributed component of UMS stars 
without an equivalent counterpart of PMS stars. 
This diffuse and "UMS dominated" component is 
clearly inconsistent with a Myr old population and 
corresponds to a (pervasive) presence of SMC field 
stars. 



In the previous section such field population has 
been considered only as a mere intruder of the 
young sub-clusters. Nevertheless, field stars also 
retain valuable information on the average SFH 
in the region. Such history is not represented 
by either BS 90 or the sub-clusters, since these 
structures are distinct snapshots in space and time 
(Gyrs old the former, Myrs old the latter) of the 
overall history in the SMC. 

We present here a quantitative analysis of the 
field SFH. To this aim, field stars (red dots in Fig. 
[T2"j) have been chosen as isolated as possible to 
avoid contamination from either BS 90 or the sub- 
clusters. However, given the radial profile of BS 90 
(see S07), the lower main sequence (fainter than 
V « 22), as well as the red giant branch (RGB) 
and the red clump (RC) may still suffer of a resid- 
ual contamination, potentially leading to an over- 
estimate of the field star formation between 4 and 
5 Gyr. 

The resulting CMD, containing about 16,000 
stars, is shown in the top panel of Fig. [131 In 
order to quantify the corresponding SFH we com- 
pared the observed CMD with an extensive set 
of simulated histories. In such simulations, the 
model behavior is changed according to the set 
of initial parameters, namely the star formation 
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law and rate SFR(t), the chemical enrichment law 
Z(t), the initial mass function IMF, the fraction of 
binary stars, the reddening and distance modulus. 
The comparison between the observed CMD and 
the model CMDs is done through minimization of 
the Poissonian \ 2 m the number of stars in a set 
of regions (0.1 x 0.1 mag large) of the CMD (cfr. 
Cignoni fc Tosill2010t ). In addition, follow i ng th e 
approach described e.g. bv iGreggio et al. ( 1998b . 
particular attention has been paid to reproduce 
the most reliable and evident evolutionary phases, 
such as the blue loops, the subgiant and the red gi- 
ant branches, the main sequence Turn-Offs. To re- 
duce computational time the SFR is parametrized 
as a linear combination of fuzzy isochrones with 
variable duration (partial CMDs). The final un- 
certainties on the SF H are obtained with a boot- 
strap technique (see ICignoni et al 1 120091 for de- 
tails) . 

A morphological comparison between observed 
and synthetic CMDs allows to reduce the param- 
eter space. We simulated mo dels with the fol - 
lowing metallicity Z(t) (see e.g. iNoel et aL 20091) : 
Z = 0.004 for stars younger than 2 Gyr, Z = 0.002 
for stars with age between 2 and 5 Gyr, Z — 0.001 
for stars older than 5 Gyr. The adopted IMF has 
Salpeter's exponent. A 30% fraction of binaries is 
assumed. Finally, our synthetic population is cor- 
rected for a distance modulus (m — M) = 18.9 
and a galactic reddening E(B — V) = 0.08. 

The first noticeable result is related to the ob- 
served UMS morphology, whose large spread in 
color is not accounted for by our models either 
with age or binaries or photometric errors. Even 
if a differential reddening of about 0.1 mag was 
effective to solve this discrepancy, the main se- 
quence blue edge in the range 21 < V < 23 would 
be still redder than our models. We suggest that 
the youngest populations in the field suffer both 
from differential reddening (which varies by about 
0.1 mag) and from a distance modulus spread of 
about 0.1 mag. 

Concerning intermediate to old generations (> 
2 Gyr), the compactness of the red clump and 
the thinness of the red giant branch rule out any 
hypothesis of differential reddening. The magni- 
tude dispersion of the sub-giant branch is consis- 
tent with a distance spread of at least 0.2 mag 
((m — M) = 18.8 — 19.00), in good agreement 
with results obtained for intermediate age star 
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Fig. 13.— Top panel: CMD for the selected field 
stars. Bottom panel: Best synthetic CMD. 



clusters ( Glatt et al. 20081 ). A suggestive scenario 
may be summarized as follows: 1) the "old" field 
belongs to an extended halo, evenly distributed 
around the average distance (m — M) = 18.9 and 
marginally affected by differential reddening; 2) 
the "young" field is localized on the farther side of 
the SMC and suffers differential reddening. 

Using these assumptions on distance and red- 
dening, we have proceeded to recover the best 
combinations of partial CMDs leading to the mini- 
mum x 2 ■ The resulting CMD and the correspond- 
ing SFH are shown in the bottom panel of Figure 
[T3]and in Figure [T4l respectively. The larger errors 
in the recent activity are mainly due to differential 
reddening. In terms of star-counts: 1) upper and 
the lower main sequence (down to V ~ 25) are 
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well reproduced; 2) the number of red clump stars 
is always slightly over-predicted, by about 30%; 
3) Blue loop stars are always under-predicted, by 
- 20 - 30%. 

The field SFR is increasing from 13 Gyr ago 
up to now (see Fig. [T4|) . We estimate the av- 
erage rate density in the last 100 Myr to be 
1.4 x 10~ 8 M yr" 1 pc~ 2 (2.5 x 10" 8 Mq yr' 1 p C - 2 
when extrapolated using a Salpeter IMF down 
to 0.1 M©) and the average rate density over the 
whole 13 Gyr lifetime to be 1.5xl0 -9 Mq yr~ x pc~ 2 
(2.7 x 10~ 9 Mq yr -1 pc~ 2 when extrapolated). 
However, about 60% (38000 Mq) of the stellar 
mass was assembled in the earliest 8 Gyr, with 
36% (21000 M ) formed between 7 and 5 Gyr ago. 
For comparison with another SMC star forming 
region, NGC 602, in the period between 3-13 Gyr 
ago the field around NGC 346 produced at least 
six times more mass per pc 2 th an the field around 
NGC 602 (jCignoni et alj|2009t ). Finally, it is also 
comforting to note that no significant SFR peak 
in the more recent period 5-4 Gyr ago is found, 
showing that the contamination from BS 90 is 
actually minimal. 

8. Discussion and Conclusions 

The analysis of the CMDs and of the spatial 
distribution of the stellar populations in the region 
of NGC 346 leads to interesting results, both for 
the clustered and the diffuse components. 

Diffuse component 

For the diffuse component we find that the ma- 
jority (60%) of the mass astration occurred ear- 
lier than 5 Gyr ago, with a peak between 5 and 
7 Gyr a go. This is at var i ance w ith the sugges- 
tion by lHarris fc Zaritskvl ( 2004 ) of a long pe- 
riod of almost no activity from 3 to 8.4 Gyrs ago, 
but in agreement with the SFH of other SMC re- 



Noel et a l. 2009, Sabbi ct al. 2009). We recall that 



gions derived by other authors (e.g. [ Dolphin et al 
20011 iMcCumber et alj|2005l ICignoni et al.l l2009. 



our photometry and those of all the latter authors 
reach the oldest MS turn-off and therefore a look- 
back time of 13 Gyr, while the shallower data by 
Harris & Zaritsky do not allow to properly de- 
scribe epochs earlier than a few Gyr. 

Over the SMC lifetime the average rate den- 



sity in the field surrounding NGC 346 is 1.5 x 
10~ 9 A/© yr^ 1 pc~ 2 . This value is in agreement 
with those derived bv lNoel et al. (2009) for several 
SMC regions. Its SFH is qualitatively similar to 
that around the other star forming r egion we stud 



ied w ith HST/ACS data, NGC 602 (ICignoni et al 
2009h . However, the activity in the NGC 346 re- 



gion is significantly higher (by a factor of 6) than 
in NGC 602. These results are not surprising given 
the much closer proximity of NGC 346 to the SMC 
center. The rate of star formation in the NGC 346 
region (and in NGC 602 as well) appears to have 
increased in the last tens of Myr from a relatively 
lower and steady regime. We must remind, how- 
ever, that the SFH recovered at very old epochs is 
averaged over timescales of Gyr and short episodes 
of higher activity may remain undetected. 

In the last 100 Myr the average star forma- 
tion rate density in the field is about 1.4 x 
10~ 8 Mq yr^ 1 pc~ 2 , a factor of 2-10 higher than i n 
the solar neighborhood (cfr. iTimmes et al. 19951) . 
two orders of magnitude higher than in nearby 
late- type dwarfs, and similar t o the quietest cases 
of Blue Compact Dwarfs (see iTolstoy et al.l 120091 
and references therein). In this respect, the SMC 
region around NGC 346 appears as currently quite 
active, but we should bear in mind that its current 
rate might be somewhat overestimated by the pos- 
sible presence of very young stars member of the 
actual star forming region itself. For comparison, 
the peak rate in SC-1 is 2 x 10 -5 Mq yr^ 1 pc~ 2 , 
three orders of magnitude higher than the field 
average value. 



Clustered component 

For the clustered component, we have found that 
the bulk of the stars in Group I sub-clusters 
are well consistent with a major star forming 
episode started about 6 Myr ago and lasted about 
3 Myr. After that, their star formation activ- 
ity has proceeded a t a low er rate. As suggested 
by IPalla fc Stahlerl (|2002l ) for the stellar group 
Taurus- Auriga, massive stars, dispersing the gas 
that would have been part of the new generations, 
may be the culprit for such "deceleration" in the 
most recent star forming activity. 

Our models do not fully explain the large color 
spread shown by PMS stars and the large observed 
number of PMS stars. The PMS shows indeed a 
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Fig. 14. — SFH of the SMC field around NGC 346. For sake of clarity the most recent 500 Myr are zoomed-in 
in the upper right onset. The uncertainty on the SFR is also shown for each age bin. 



redward broadening which is not reproduced by ei- 
ther our youngest isochrones or photometric scat- 
ter. iHennekemper et al.l (|2008l ) suggested that a 
combination of differential reddening, variability 
and binarity may account for the observed spread. 
We consider variable reddening affecting individ- 
ual PMS stars a more likely explanation (also tak- 
ing into account that our PMS evolutionary tracks 
do not include circumstellar reddening). More- 
over, our models seem to underestimate the ob- 
served PMS counts by about 40% (possibly and in 
part because our PMS tracks do not cover masses 
smaller than 0.45 Mq). We speculate that a very 
young generation of stars, so young not to have 
had time to assemble more massive stars, could 
account both for the observed excess of PMS stars 
and for the intrinsic redness of such stars; other- 
wise, the IMF must be steeper than Salpeter. 

While the origin of these discrepancies is still 
unclear, the comparison with CMDs of stars in 
different locations provides clues to complete the 
puzzle. 



Low mass PMS stars dominate the star counts 
among Group II sub-clusters, showing a ratio 
PMS/UMS stars that is at least three times higher 
than in Group I. Such a population should have 
a MS counterpart, which is however not visible. 
Among all the sub-clusters in NGC 346, those in 
Group II host the reddest PMS stars. Only a 
strong and differential reddening correction could 
reconcile the Group II PMS with the SC-1 PMS, 
but it would pay the price of exacerbating the lack 
of UMS stars. Vice versa, an age effect, with the 
Group II sub-clusters experiencing now their first 
episode of star formation, may explain both the 
PMS redness and the higher PMS/UMS ratio. If 
the NE-Spur, where the Group II sub-clusters are 
located, is indeed the youngest region of NGC 346, 
a fraction of stars may have managed to retain the 
circumstellar material and appear redder. More- 
over, the youth of the system and the peculiar 
distribution of NE-Spur stars, preferentially lo- 
cated along filaments, may determine the paucity 
of UMS. In this case the peak gas densities in the 
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filaments would be insufficient to produce massive 
stars, such as those detected in SC-1. As a con- 
sequence, the ratio of low-mass to massive stars is 
higher where the density of massive stars is lower. 

Th is trend is also noted by IPanagia et al. 
(|200d ) in the field of SN1987A. It indicates that 
star formation processes for different ranges of 
stellar masses are rather different and/or require 
different initial conditions. An interesting corol- 
lary may be that, if the youngest objects have a 
filamentary distribution (where presumably the 
gas density was higher), a critical density thresh- 
old exists below which the star formation is sup- 
pressed. The reason for this is probably related to 
the pristine conditions in the parental molecular 
clouds, like temperature or turbulence. 

As an alternative to environmental effects, 
the lack of UMS stars may reflect the build up 
timescales of different masses. More massive stars 
may form later than low-mass stars because they 
need more time to collect enough material to start 
the formation. In this case the Group II sub- 
clusters are simply too young to have produced 
massive stars. 



The analysis o f 



Gouliermis et al 



Hcn nekemper et al.l (|2008) and 
2008j) provide an independent 



support to the age hypothes is. Fitting Seiss 

isochrones, iHennekemper et al. (2008) derived for 
the NE-Spur sub-clusters ages young er than 2.5 



Myr, while iGouliermis et al.l (|2008l ) concluded 



that the star formati on there was recently t rig- 
gered (see figure 1 in IGouliermis et al.l 120081) by 
the massive proge nitor of SNR B005- 7724. We no- 
tice however that iNaze et al.l (J2002j) argued that 
the SNR should be located in front of the lumi- 
nous blue variable HD-5980 and not really within 
the NGC 346 region. 

From the point of view of the gas, the NE- 
Spur is a natural place to find new stellar gen- 
erations. Spectroscopic observations of CO prove 
the existence of reservoirs of cold molecular gas 
in the spatial regi on around SC-14 and SC-15 (see 
Rubio et al.ll2000l ) while the distribution of the ra- 
diation field at 160 nm is clearly confined to the 
Bar, and it is not detected in the Spur (except 
in SC-13). In other words, the extreme NE-Spur 
contains residual reservoirs of molecular gas that 
may be fueling the star formation. 

The Group III sub-clusters exhibit: 1) the same 



anomalous ratio PMS/UMS stars found in Group 
II (hence at variance with Group I); 2) a color 
spread and redness among LPMS that is perfectly 
in line with what is found in Group I (hence at 
variance with Group II); 3) a color spread among 
UMS that is at variance both with Group I and 
with Group II members. 

This suggests a picture where the Group III 
sub-clusters are in an intermediate state between 
those in Group II and those in Group I. If the un- 
derlying parameter is the age, these sub-clusters 
are sufficiently old that low mass PMS stars 
have already completed their accretion phase 
(hence their color dispersion is consistent with 
the "evolved" cluster SC-1), but at the same time 
so young that intermediate and massive stars are 
still approaching the UMS. As for the NE-Spur, 
another key of interpretation is the available gas 
out of which the presently observable stars were 
assembled. In this picture, the gas density in the 
NW-Bar was such that intermediate and massive 
stars were formed more slowly than in the SE- 
Spur. 

These findings appear to support the view 
that low-mass stars form more "easily" than 
massive ones either because they need less gas 
density or lower rates of accretion. Observa- 
tio nal support for this wa s alrea dy presented e.g. 

by iRuppert fc Zinnecker ( 2009 ). Likewise, the 

Spitzer detection ( Simon et al. 20071 ) of proto-OB 
stars in Group II and Group III sub-clusters may 
explain the apparent lack of UMS stars in terms 
of age. 

As comprehensively reviewed by lZinnecker fc Yorkd 
(|2007l ) no consensus has been reached yet on which 
is the most likely process for massive star forma- 
tion: monolithic collapse in isolated cores, com- 
petitive accretion in a protocluster environment 
or stellar collision and mergers in very dense sys- 
tems. NGC 346 seems to favor one of the latter 
two (or both) but we should wait for further high- 
resolution observations at longer wavelengths (e.g. 
with WFC3 on board of HST, but also ALMA and 
JWST, and eventually with ELTs) to get a bet- 
ter insight in its SF process based also on its still 
embedded, not visible, youngest objects. 

In the conditions described above an assessment 
of the IMF in the NGC 346 region is risky. There is 
no doubt that in several sub-clusters the number 
of existing massive stars is definitely lower than 
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predicted by a Salpeter IMF. ISabbi et al.l (|2008l) 
already pointed out that in the region massive 
stars are underrepresented. We speculate that the 
paucity of UMS in NGC 346 has a double ori- 
gin: the maximum mass populating Group II sub- 
clusters is a consequence of the radiation feed- 
back, while the maximum mass populating the 
PMS sub-clusters is mainly a matter of youth. 
From a general point of view, our result goes in 
the same direction of a relation between the mass 
of the most-massive star in the cluster and the 
mass of its parent star c luster as suggested by 
Weidner fc Kroupa ( 20061 ). On the other hand, 
if indeed clusters form in an ordered fashion pro- 
ducing first low mass stars and then proceeding to 
assemble more massive stars till a maximum mass 
is born whose feedback halts the collapse, how is 
it possible to explain the existence of sub-clusters 
hosting only PMS stars? We expect Alma, Hcr- 
schel and JWST to provide the necessary clues. 
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